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It is argued that QCD might solve the strong CP problem on its own. To test this idea, a lattice simulation 
suggests itself. In view of the difficulty of such a calculation we have, as a first step, investigated the problem in 
the CP'^ model. The CP'^ model is in many respects similar to QCD. In this talk 1 shall present some first results 
of our calculation. Among other things it is shown that the model has a first order deconfining phase transition 
in 9 and that the critical value of 9 decreases towards zero as /? is taken to infinity. This suggests that 9 is tuned 
to zero in the continuum limit. 



1. INTRODUCTION 

The discovery of instantons in non-abelian 
gauge theories [1] has led to the solution of the 
Ua(1) problem in QCD [2]. Yet it created another 
puzzle: the strong CP problem. 

Instantons represent tunneling events between 
vacua of different winding number n, which we 
denote by \n) . They cause the perturbative vac- 
uum |0) to become quantum mechanically unsta- 
ble. The proper vacuum states, which are stable 
under any gauge invariant operation, are 

|(?) =^exp(i(?n)|n). (1) 

n 

These so-called vacua correspond to the action 
Se=S- i0Q, (2) 
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where S is the standard action and 

Q = j d'^^<^p.i'pcrT^'^Ffj_„Fpa e 2 (3) 

is the topological charge. By convention I will 
choose 9 to lie in the interval [0,27r). The ad- 
ditional term in (2) breaks parity, time-reversal 
invariance and CP. In the presence of quarks 
the 9 angle can be rotated into a phase of the 
quark mass matrix M and vice versa, due to the 
Ua(1) anomaly. The effective angle is in this case 
9 = 9 -\- arg det M . A priori 9 is a, free parame- 
ter. Experimentally one obtains the upper bound 
9 < 10-3 [3]. 

Why is 9 so small? Three possible explanations 
have been offered, all of which involve imposing 
global symmetries on the low energy physics. One 
explanation is that the u quark is massless. In 
this case the action is invariant under a chiral 
transformation, and so the 9 parameter can be 
rotated away. However, I do not think that a 
theory with a massless charged particle does ex- 
ist [4]. In a second explanation [5] an additional 
U(l) symmetry is introduced which dynamically 
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tunes 6 to zero. This implies the existence of a 
light particle, the axion [6]. If it exists, its mass 
must lie between 10~® and 10~^ eV [7]. The third 
explanation is based on introducing a discrete CP 
symmetry which is broken spontaneously. The 
parameter 6 is then calculable and can be made 
small enough to be consistent with experiment 
[8]. 

In addition, several authors have sought a solu- 
tion of the strong CP problem wtthm QCD, that 
means without the introduction of new symme- 
tries and particles [9]. It is almost certain that the 
QCD vacuum will undergo fundamental changes 
when 9 is taken different from zero [10]. Some 
time ago 't Hooft and Mandelstam have conjec- 
tured [11] that the confinement phenomenon can 
be understood in terms of a color magnetic su- 
perconductor, in which color magnetic monopoles 
condense and color electric charges, i.e. quarks 
and gluons, are confined by a dual Meissner ef- 
fect. More recently this picture of the QCD vac- 
uum has been successfully tested by lattice sim- 
ulations [12]. Not only does one find evidence 
for monopole condensation, but it could also be 
shown that the monopole currents around a fiux 
tube satisfy a dual London equation [13]. In the 9 
vacuum these monopoles acquire a color electric 
charge of the magnitude 6'/27r [14]. For 9^0 one 
would then expect [15] that the long-range color 
electric forces are screened by monopoles and that 
confinement is lost. Very likely this will cause the 
breakdown of asymptotic freedom as well. 

Probably this is not the only way of under- 
standing the QCD vacuum. For example, it might 
also be possible to describe the infra-red proper- 
ties of the QCD vacuum in terms of instantons. 
But I would expect that one will find the same 
dependence on 9. Indeed, in a model of the QCD 
vacuum based on semi-classical ideas it has been 
argued that in the 9 vacuum the long-range color 
fields are screened by instantons [16], thus lead- 
ing to similar conclusions. (A further example is 
given by the CP^ model which I will discuss in de- 
tail below. This model possess instantons but no 
monopoles, and the vacuum can be understood 
as an instanton liquid.) 



How would this solve the strong CP problem? 
If it proves to be true that confinement is lost for 
9^0, then 9 = is the only choice. The problem 
then is to understand why the confinement mode 
is chosen by nature. The answer could be that 
the theory admits a continuum limit only for 9 = 
0. 

To test this possibility, the natural way to pro- 
ceed is to simulate the theory on the lattice. 
What makes this calculation very difficult, how- 
ever, is the fact that the action is complex for 
non-vanishing values of 9. As a result, standard 
lattice techniques are not immediately applica- 
ble. This has led us to investigate the problem in 
a simpler model first. 



2. THE CP'^-i MODEL: A 2-D IMAGE 
OF QCD 

A model, which in many respects is similar to 
QCD, is the CP^~^ model in two space-time di- 
mensions. The common properties include (i) the 
existence of instantons and a vacuum angle 9, (ii) 
asymptotic freedom, (iii) a dynamically generated 
mass gap, (iv) dimensional transmutation and 
(v) confinement. A linear potential is of course 
not difficult to achieve in two dimensions, but in 
the CP^~^ model confinement arises without the 
presence of any fundamental gauge fields which is 
remarkable. 

The CP^~^ model describes N-component, 
complex scalar fields Za(x) of unit length: 

Za(x)Za(x) = I, a = I,- ■ ■ ,1^. (4) 

Out of these fields one may construct composite 
vector fields 

= ^^l Za{x). (5) 

Then the action [17] can be cast into the intuitive 
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form [18] 



S = (3 I d xDf,zJx)Df,Za(x), 



(6) 



where = 5^ + iA^. It follows that (6) is in- 
variant under the gauge transformation z'^(x) = 
exTp(ia(x))za(x). Thus the model corresponds to 
a set of charged scalar fields interacting minimally 
with a composite gauge field. The topological 
charge is given by 

£xie^^D^z^{x)D^Za{x) (7) 



Q 



27r 
1 

2^ 
1 

2^ 



d^3:€fj_i,dfiAi,{x) 

(PxFqi. 



(8) 
(9) 



This means that a field configuration with a non- 
trivial topological charge can also be seen as a 
configuration with a background electric field. 

Confinement is directly related to the 9 depen- 
dence of the theory, as we shall see in a moment. 
The 9 dependence is essentially described by the 
partition function 



z{9) = Y,^.Mm)piQ), 

Q 



(10) 



where 



p{Q) 



j[VzVl]Q8{\z\^-l) exp(-g) 
]vzVz8{\z\'^ - 1) exp(-5) 



(11) 



The subscript Q means that the path integration 
is restricted to the charge sector indicated. The 
associated free energy per space-time volume V 
is 



F{9) 



\nZ{9). 



(12) 



The connected moments of the topological charge 
distribution are obtained from F(9): 



{o\Q"\o)c 



,d"F(9) 



(13) 



In particular, the average charge density is given 

by 
1 



V 



{9\Q\9) 



-W) 

.dF{9) 



(14) 
(15) 



Consider now a pair of static external particles of 
charge e and — e in units of the intrinsic charge. 
The potential between the two particles is derived 
from the Wilson loop formula 



W{A,9) = {9\eM^e i dx^A^)\9), (16) 

Jc(A) 

where A is the area enclosed by the contour C{A). 
From (8) and (15) it then follows that 

W(A,9) = exp(-A / d9'—^) (17) 



= exp(-A[F(9 + 2TTe) - F(9)]) (18) 

which for 9 = and small e implies confinement 
by a linear potential. However, when e = 1, i.e. 
when e is equal to the intrinsic charge, the long- 
range force is absent because F(9) is periodic in 
9: F{9) = F(9 + 2tt) [19]. The physical reason 
for that is that the external charges are screened 
by pair-produced z particles. This is exactly the 
same situation as in QCD. The string tension is 
given by 



a(e,9) = F(9 + 2TTe) - F(9). 



(19) 



Note that the string tension is not invariant under 
charge conjugation, i.e. a(e,9) ^ a(—e,9), for 
9 ^ 0,17 due to the fact that the vacuum is an 
eigenstate of CP only for 9 = 0,7r. Thus, if a 
particle of charge e is to the left of an anti-particle 
of charge — e it will feel a different force than if it 
is to the right of the anti-particle [20]. 



3. LATTICE SIMULATION WITH 0^0 

The CP^~^ model has been investigated by 
large-N [18] and instanton methods [21] in the 
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past. In the former approach no spectacular 6 de- 
pendence was observed in leading approximation. 
But this is not surprising. The higher moments 
(n > 2) of the topological charge distribution (13) 
are zero in this approximation which eliminates 
quantum effects to a great extent. 

We have chosen to examine the CP^ model in 
closer detail [22]. On the lattice the action takes 
the form 

S = -2(3j2\M^)^a(x + i:i)\'' (20) 



are selected randomly [25]. This is supplemented 
by an additional acceptance criterion: if the new 
charge is in the same set the configuration is ac- 
cepted and the new charge recorded; if the new 
charge is outside the set the change is rejected 
and the old charge recorded. Furthermore, a trial 
charge distribution which is approximately equal 
to the true distribution is incorporated. 

The probability function p(Q) is a steeply 
falling function of Q. With this method we were 
able to compute p(Q) over a range of more than 
twenty orders of magnitude. 



(a = l,---,4). The boundary conditions are 
taken to be periodic. In accordance with (5) one 
defines link fields 



Za(x)Za(x + fl) 

\zi,(x)zi,(x + fl)\ 
exp(iAi_,(x)). 



(21) 
(22) 



The topological charge is derived from (21) by 
forming parallel transporters around plaquettes, 



Un = exp(i_Foi), -t < Fqi < tt, 



and writing 



(23) 



Q 



(24) 



It is easy to see that Q is an integer. For this 
definition of Q it is expected that the model does 
not give rise to dislocations [23]. Dislocations are 
lattice artifacts on the scale of the order of the 
lattice spacing which cause the free energy F(6) 
to diverge in the continuum limit. 

We compute Z(6) by computing the probabil- 
ity function p(Q). This is done by the following 
method [24]. We divide the phase space into over- 
lapping sets of five consecutive charges. In each 
of these sets the z fields are updated by a com- 
bination of Metropolis and overrelaxation steps, 
z' = exp(i(/)aAa)z, where the ?7(4) generators Xa 



4. FIRST PROMISES 

We are ready now for a quantitative test. If 
our idea is correct, we should find a first order 
phase transition in 6 from a confining phase to a 
Higgs or Coulomb phase. On a finite lattice and 
at a finite value of /3 the phase transition would 
occur at a value 6 = Oc(p, V) > 0, where V now 
is the lattice volume. Only on an infinite lattice 
and at /3 = cxd, i.e. in the continuum limit, can 
we expect that 9c = 0. 

A first order phase transition will manifest it- 
self in a kink in the free energy, F(6), as well as 
in a discontinuity in the first derivative of the free 
energy, dF(9)/d9. In Fig. 1 I show F(9) on the 
V = 64^ lattice at /3 = 2.7. ^ The correlation 
length for this value of /3 is ^ 8.8 [26]. For com- 
parison I also show the prediction of the large-N 
expansion to leading order [18], which turns out 
to be F(9) = c9^, where the constant c has been 
fitted to the lattice data at small 9. ^ We see a 
distinctly marked kink at 9 = 9c ^ 0.5 tt: while 
F(9) increases roughly proportional to 9^ up to 
9 = 9c, we find that F(9) is constant (within 
the error bars) for 9 > 9c. The latter result is 
rather remarkable, as it means that all deriva- 
tives of F(9) vanish for 9 > 9c. In other words. 



^The reader who recalls Ref. [22] might notice a slight 
change in the shape of F{0). The explanation is that in 
the meantime we have increased our statistics by an order 
of magnitude. 

3This gives c= {l/2)d^F{e)/de^\e=o = (1/2) (OlQ^ |o)/y. 



4 




0.0 0.2 O.M 0.6 0.8 1.0 

e/7T 



Figure 1. The free energy F(6) as a function of 
e on the V = 64^ lattice at /3 = 2.7. The solid 
curve is the prediction of the large-N expansion to 
leading order. Only the first half of the 6 interval 
is displayed. In the second half of the interval 
F(0) = F(2tt-0). 



the theory becomes independent of 6. The large- 
N expansion, on the other hand, predicts a phase 
transition at 6 = it. In Fig. 2 I show the first 
derivative of the free energy. According to (15) 
this quantity is identical to the topological charge 
density, q(0), which again can be interpreted as a 
background electric field (cf. (9) and (24)). We 
see that q(6) increases almost linearly with 6 up 
to 6 = 6c, where it jumps to zero and then stays 
zero over the rest of the interval (again within the 
error bars). Thus the phase transition is accom- 
panied by a collapse of the background electric 
field which hints at a transition from a confining 
phase to a Higgs phase. 

Before we can make any statements about the 
continuum, we have to extrapolate the lattice re- 
sults to V = CO and /3 = cxd. For a first order 
phase transition we expect 

9rXP,V)-9rXP,^)ocV-^ (25) 
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Figure 2. The topological charge density q(6) as 
a function of on the V = 64^ lattice at /3 = 2.7. 
Only the first half of the 6 interval is displayed. In 
the second half of the interval q(6) = —q(2'iT — 9). 



for fixed [3. In Fig. 3 I show the results for 
OciP, V) for a variety of lattice volumes ranging 
from V = 72^ to V = 28^ and for two values of 
P, P = 2.5 and /3 = 2.7. For /3 = 2.5 the cor- 
relation length is ^ 4.5 [26]. I have chosen to 
plot 6'c(/3, V) as a function of . We see that 
for both values of [3 our data fall on a straight 
line in accordance with a first order phase transi- 
tion, thus confirming our earlier statement about 
the order of the transition. This allows us to ex- 
trapolate the lattice results to the infinite vol- 
ume. We obtain ejiT = 0.32(2) for /3 = 2.5 and 
Bc/iT = 0.18(3) for l3 = 2.7. These values are 
to be contrasted with the analytic result in the 
strong coupling region 9c = it [27]. 

Three values of /3 are not enough to make a 
precise extrapolation to /3 = cxd. But the trend 
is evident: 6c depends strongly on /3, and 6c de- 
creases towards zero as we approach the contin- 
uum limit. Our resuls are consistent with a decay 
like 6c oc 1/^. This suggests a phase diagram of 
the form sketched in Fig. 4. The horizontal line 
at 6 = IT is the prediction of the strong coupling 
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Figure 3. The critical value Oc(p, V) as a function 
of for two values of /3. The lines are a linear 
fit to the data points. 



calculation. It is a line of first order phase tran- 
sitions on which CP invariance is spontaneously 
broken. Whether 6c is indeed exactly zero at 
/3 = cxD we cannot say at the moment. But our 
results so far are promising. In the near future 
it should be possible to repeat the calculation for 
a correlation length that is up to a factor of five 
times larger than that at /3 = 2.7 [25, 26]. 

To show that the phase transition is indeed a 
deconfining phase transition, let me discuss the 
string tension now. In (19) I have derived a for- 
mula that relates the string tension to the free 
energy. According to this formula it is sufficient 
to know F(6) over a wide enough range of 9. This 
we do already. To keep the discussion simple, 
I shall restrict myself to the case of small frac- 
tional charge e. This matters because the vac- 
uum will change its properties if it is exposed to 
a strong external electric field, as we have seen 
before. The 6 dependence of the string tension 
can then be read off from Fig. 1. From that fig- 
ure and (19) it follows that at 6 = the string 
tension is a(e,0) k, c(27re)^ = 2i:'^e^xt, where 

= (0|(5^|0)/y is the topological susceptibility. 



Figure 4. The phase diagram over the full interval 
Q < 9 < 27r. The horizontal line at = tt is the 
strong coupling prediction, and the symbols are 
the lattice data extrapolated to infinite volume. 

This relation becomes exact in the limit e ^ 0. ■* 
If we now increase 9, the string tension will in- 
crease until 9 reaches 9^ — 27re, where it will start 
to decrease again. Finally, at 9 = 9^ the string 
tension will be zero. And provided the external 
charge is small enough, it will remain zero up to 
9 = IT. This result is a consequence of the prop- 
erty that F(9) is constant for 9 > 9^. Let us now 
look at the data for the string tension directly. 
In Fig. 5 I show <T{e,9) for two charges, e = 0.1 
and e = 0.2, on the V = 64^ lattice at /3 = 2.7. 
Remember that on this lattice 9c was 0.5 tt. 
We clearly see that the string tension vanishes 
for 9 > 9c (within the error bars). 

We have also investigated the /3 function and 
its dependence on 9. We find evidence that the 
l3 function vanishes at the phase transition, indi- 
cating that asymptotic freedom might also be lost 
in the so-called Higgs phase. Because of lack of 
space I cannot go into detail here. A full account 
of our work will be given elsewhere [28]. 

^Note that this is precisely the resuU of the large-N ex- 
pansion to leading order [18]. 
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e=0.2 



e=0.1 



that the theory is not asymptotically free in the 
Higgs phase. I would therefore not expect to find 
any other ultra-violet stable fixed point. 

So far we have focussed all our efforts on the 
CP^ model. The algorithm we have employed has 
proven to be very efficient. I am confident that 
one will be able to solve the strong CP problem 
one day also in QCD, though the demand in com- 
puting power will be vast. As a first step in that 
direction we are currently investigating the prob- 
lem in the SU(2) Yang-Mills theory. 



e/7T 
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Figure 5. The string tension <T{e, 9) as a function 
of 6* for e = 0.1 and e = 0.2 on the V = 64^ lattice 
at l3 = 2.7. Only the first half of the 6 interval 
is displayed. In the second half of the interval 
a(e,e) = a(-e,2TT - 6). 



5. SUMMARY AND OUTLOOK 

Let me now summarize our results. We have 
found a phase transition in 9 from a confining 
phase to a non-confining phase, presumably a 
Higgs phase. This by itself is already a remark- 
able result which has interesting applications. 
The critical value of 9 turns out to vary strongly 
with [3. While 9^ = i: in the strong coupling re- 
gion, 9c falls off rapidly in the direction of zero 
as [3 is taken to infinity. This behavior is summa- 
rized in the phase diagram shown in Fig. 4. The 
important question that remains to be answered 
is whether the line of first order phase transitions 
will indeed end in the point 9) = (cxd, 0). If 
this proves to be the case, then 9 = Q is the only 
point at which the continuum limit can be taken, 
at least in the confining phase. This would resolve 
the strong CP problem. As far as other possible 
fixed points are concerned, there are indications 

consider only the lower branch of the phase transition 
line here. 
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